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Abstract 

We discuss a class of transfer matrix built by a particular combination of iso- 
morphic and non-isomorphic GL(N) invariant vertex operators. We construct a con- 
formally invariant magnet constituted of an alternating mixture of GL(N) "spins" 
operators at different order of representation. The corresponding central charge is 
calculated by analysing the low temperature behaviour of the associated free energy. 
We also comment on possible extensions of our results for more general classes of 
mixed systems. 
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1 Introduction 



One of the most useful methods of constructing an integrable one dimensional quantum 
spin chain has been to find solutions of the Yang-Baxter fl], equations. In the context 
of the lattice models, such solutions define vertex operators acting on the tensor product 
of two vector spaces V ® h a of the local states on the horizontal and vertical lines of a 
given site a of the two dimensional lattice. The complete Hilbert space of the model on a 
lattice of L sites is Yia=i ®h a . The vector space V is an auxiliary space which is useful in 
the definition of the associated transfer matrix. Defining R^ a (/j,) as such vertex operator, 
the corresponding transfer matrix can be expressed by [?, || 

T(»)=Tr v [RlMRl L _M---Rl 1 ( f i)} (1) 

where the trace is carried out in the auxiliary space V and [i is a variable which parametrizes 
the Yang-Baxter solution. 

A class of important solutions of the Yang-Baxter equation is that isomorphic on the 
horizontal and vertical spaces (V = h a ) and rational on the spectral parameter fi. One ex- 
ample is the vertex operator composed of generators invariant by some of the semi-simple 
A, D, E Lie algebra at certain order k of its representation. The main feature of these 
solutions is that they are believed to define conformally invariant quantum spin chains 
[|], |28j which realize the class of universality of Wess-Zumino-Witten-Novikov (WZWN) 



U field theories with topological charge k. Recently, de Vegaand Woynarovich M have 
pointed out that other interesting classes of conformally invariant models can be con- 
structed. These theories are obtained by combining isomorphic and non-isomorphic vertex 
operators invariant by some group symmetry at different order of representation. As a 
typical example one can consider an alternating combination of isomorphic vertex at even 
sites and of non-isomorphic operator at odd sites. Considering that the isomorphic (non- 
isomorphic) vertex operator acts on the vector space ® hffl (l/( fc ) <g> )) of order 
k and k\ the composed theory described above will lead to a mixed magnet chain with 
"spins" operators of order k and k' . Following this approach, de Vega and Woynorovich 



U have constructed an alternating anisotropic Heisenberg chain of spins 1/2 and 1 and 
analysed its behaviour in the thermodynamic limit. More recently, the central charge of 
the conformally invariant isotropic 577(2) mixed chain has been computed independently 
in refs. [0, |J for spins 1/2-1 and 1/2-S, respectively. 

Its seems interesting to investigate the critical properties of conformally invariant mixed 
spins chain for a more general class of group symmetry and its possible relations with 
WZWN field theories. In this paper we study an isotropic alternating GL(N) model in 
which the isomorphic vertex is in the fundamental representation (k = 1) and the non- 
isomorphic operator is defined in the tensor space product of the fundamental and the 
order k of representations. By using the thermodynamic Bethe ansatz, we have calculated 
the associated central charge in the case of a conformally invariant mixed GL(N) system. 
It is noted that this central charge can be decomposed in terms of the conformal anomaly 
of two WZWN theories with different topological charges k = 1 and k = k — 1, respectively. 
We also discuss the generalization of our results in the case of semi-simple Lie algebras at 
arbitrary symmetric representation. 

This paper is organized as follows. In sect. 2 we define the GL(N) mixed magnetand we 
discuss its diagonalization by the quantum inverse scattering (QIS) formalism. In sect. 3 
we use the thermodynamic Bethe ansatz (TBA) approach in order to investigate the low 
temperature behaviour of a conformally invariant mixed GL(N) system. Sect. 4 is devoted 
to our discussion on possible extensions of the results of sect. 3 to other Lie algebras. Sect. 5 
contains our conclusions. In Appendices A and B we summarize some details of the QIS 
approach and we present extra numerical checks of the finite size behaviour for the ground 
state energy, respectively. 

2 The mixed GL(N) integrable model 

A rational GL(N) invariant vertex operator has been found by Kulish, Reshetikhin and 
Sklyanin || on their study of group invariant solutions of the Yang-Baxter relation. The 
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GL(N) non-isomorphic vertex defined on the fundamental and on the symmetric order k 
of representation is given by the expression [[J 

RlM =/*-(*- 1)/2 + E ^' ® ^ (2) 

where e 4J are the generators of GL(N) in the fundamental (order k) representation. 

For the fundamental representation the matrix elements of e %3 are (e y ')jy = SaSjk- We also 
notice the identity ij[ a (0) = V, where V is the permutation operator VV ®h a = h a ®V. 

The GL(N) mixed system is defined in terms of its transfer matrix of alternating 
isomorphic {k = 1) and non-isomorphic vertex operators by 

Ti,k(n) = Tr vW [r lik (fi)}, n tk (fj) = R\MRk,L-x(^) • • • R\,MR\,M (3) 

where the matrices product and the trace are defined in the auxiliary space V {1) of the ele- 
ments e %3 and Tx^ip) is the so-called monodromy matrix. The associated one dimensional 
quantum Hamiltonian acting on the Hilbert space Y\a=i ^« * s determined as the logarithm 
derivative of the transfer matrix, Hx,h — [ln(T l fc (/j,)] 1^=0- By using this relation and 
some properties of the GL(N) group we find 

N N 



J L_1 
H hk = g E 

n=odd 



E 4-i{Et, E^e-U - (k - 1) E(eti^4V! + eUEfrti) 

i,j,l,m 



\2 N N 



+ { -^r L E 4-i<&i + Y^LM + e^: +1 ) -(k-r 



h3 hO 



— (1 + T^-) (4) 
2 v k+r v ; 



where J = 4 J / {k + l) 2 and we have conveniently added an extra constant — J\{\ + j^zi)- 
In this paper we are interested in the antiferromagnetic properties of this model, and 
therefore we have chosen J = 1. Setting k = 1 in Eq.(£|) we reproduce the results of 
refs.|f[3|, pT0|] and for N = 2 and k = 2S we recover previous calculations for the mixed 



Heisenberg chain || |], |7|] . 

The diagonalization of the transfer matrix Tx^in) (or the one dimensional Hamiltonian 
Hx,k) follows from the generalization of the quantum inverse scattering (QIS) approach || 



applied to multi state vertex models [10, O, U2fl. In the QIS construction the definition of 



the pseudo vacuum and the block form of the monodromy matrix play as two important 
ingredients of this method. In our case the pseudo vacuum |0 > is defined by 

|0 >= |0 >\ ®|0 >2 ® ■ ■ ■ ® |0 >£_i <8>|0 >i (5) 

where |0 >Jj is the vector of highest weight of the GL(N) algebra at order k of the 
symmetric representation acting on the site a, namely E%>\0 > k a = fc^i^ilO An 
important property of this state is that Tijfe(//)|0 > has the following block triangular form 

T ljh (fi) >= |0 > (6) 

V 5ij(fi) L/2 (fi — (k — 1)/2) L/2 J 

where = n^) 1 ^ 1 , i = 1, • • • , N - 1. 

Following the QIS machinery jlO], 11, [12| a certain linear combination of the states 



B^d^l) • • • JB 3 ' n (//^)|0 > can be considered as a basis for the eigenstates of the transfer 
matrix, provided that the parameters {fi\, ■ ■ ■ , jj^} satisfy a set of non-linear equations 
denominated Bethe ansatz equations. The technical steps are fairly parallel with those of 



refs. [10, 11, 12] and in appendix A we have collected some of the details. Defining the 



convenient shift fi r j = iX'j — r/2, the Bethe ansatz equations are given by 

( x)-i/2 \ L/2 ( X]-ik/2 \ L/2 = _ Ml Xj-Xj-i ft A) - Aj + i/2 
\\)+i/2) \X) + ik/2) ~" l\ X) - A 1 + % l\ X) - A 2 - i/2 

M r \r \r ■ M r+1 \r+l \r+l ■ in Af" 1 \r-l \r-l -In 
Aj - A l - 1 = Aj — A l - 1/2 Aj — A[ - 1/2 

t\ ^ ~ A r + i l\ A f 1 " Ar+1 + L\ X T l - Ar_1 + z /2 ' 
where r = 2, • • • , N — 1 and we define M N = 0. The eigenvalues of the Hamiltonian H%,k 
are parametrized by 

A/ 1 i 



U (A}) 2 + (1/2) 2 

In principle, the same construction discussed above can be carried out for the ' 'dual" 
one dimensional Hk,i spin chain. The main task is to find an isomorphic GL(N) invariant 
operator in the symmetric representation of order k. Using the approach of ref. 
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H, Johannesson ]14| has explicitly exhibited such operator. Hence, analogously to Eq.fl3|), 



the associated transfer matrix T k ^(fi) is expressed by 

T M ( M ) = Tr v(k) [r k M], r M (//) = ^(^t-iM ■ ■ ■ R\MR k k M (9) 

where is the space of the matrices of GL(N) at order k of representation. 

An important property of Ti.fc(/i) and T^i(fi) is their commutativity for arbitrary 
values of the parameters \x and //. This property follows from the "mixed" Yang-Baxter 
relation satisfyed by the non-local vertices^ R\{p) and R\{p). As a consequence, the 
Hamiltonians H± tk and H kj \ can be simultaneously diagonalized and their eigenspectrum 
are parametrized by the same Bethe equations, i.e., Eqs. (]?]). However, the eigenenergies 
of H ky i [14|] are expressed in terms of the variables {A]} by a different function of that of 
Eq.(||), namely 

m 1 k 

Bfc ' 1 = -g(A]FT(W (10) 

At this point, it is important to remark that the transfer matrix Tx tk (fi) and its "dual" 
T k} i(n) are not rotational invariant in the horizontal/vertical space of states. A simple way 
of defining |J a "mixed" symmetric transfer matrix T sym (n), preserving the rotational 
invariance, is by formally multiplying these two transfer matrices 

T^/i) = T lifc (//)T M (//) (11) 

Due to the commutativity between Ti k (fi) and T k i(fi) the spectrum of T sym ({i) is 
clearly parametrized by the same Bethe ansatz equation and the eigenvalues of the corre- 
sponding one dimensional Hamiltonian H sym are added, 

M 1 -, M 1 u 
E sym \- L _ V - (12) 

(A}) 2 + (1/2) 2 £ (A*) a + (*/2) a 1 ] 

In this sense Eqs.(|7], || [l(], |12| ) define three families of possible spectrums parametrized 
by a single Bethe ansatz equation. For instance, in the case of small size L, one only needs 



1 The vertex R\ (jj) acts in the non local space V^ k ) <g> h. 
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to solve Eq. (|7|) and compare it with the exact solution of the spectrum of in order to 
investigate the structure of the variables {A^} which are going to parametrize the whole 
spectrum of all these three families of models. In general, we remark that a certain solution 
{X r j} of Eq.([7D will not necessarily produce the same i th state ordered in energy for all these 
models. The ground state, however, is characterized by the same structure of {A^} for all 
these three families. In the thermodynamic limit, L —>■ oo, the ground state solution {A^} 
is composed by a mixture of real numbers rf- and a set of complex roots X r j a . The complex 
structure X^ a cluster in the so-called k-string form 



A 



J. a 



(k + l-2a), 



a: 



1,2,.-.,/: 



(13) 



where £J is a real parameter denominated center of the k-string. 

In order to calculate the ground state energy in the L — > oo limit, one has to solve 
the Bethe ansatz equations ([?]) for the variables rf- and xij. Taking into account Eq. ([T3|) , 
and after some standard manipulations we are able to compute the following values for 
the ground state energy per particle 

,k-l + 2N, ,,k+\ 



1 

'iV 

1 

'n 



V>(i)-^(i/;v) + v( 

^(l)-iP{k/N)+i)( 



2N ' 

k- 1 + 27V\ 
2N ' 



2N ' 

k + 1, 
2N ' 



k-l AT 

i=i J 



(14) 
(15) 



where if)(x) is the Euler psi function. From Eq.(12) it is clear that = + e 



k,l 
oo • 



Let us now concentrate our attention on the low-lying excitations of the symmetric 
mixed model. Since the transfer-matrix T sym (n) is rotational invariant, its associated 
quantum spin Hamiltonian is a strong candidate to be conformally invariant. We recall 
that the analysis of the critical properties in spin chains depends on the behaviour of its 
dispersion relation for low momenta p. The computation of the dispersion relation follows 
the standard formalism of perturbing the ground state structure by holes and string of 
arbitrary length (see e.g. fl^fl). The only subtle fact is that this alternating mixed system 
the total momenta is half of those considered in homogeneous models (k = 1) ||. Following 
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refs. [13| , |l4 > || we find that the dispersion relation for all branches of excitations |T3|, [14 



possess the same linear behaviour for the total low momenta p 

47T 

e(p) = j^P ( 16 ) 

From Eq.(|l6D the sound velocity is v s = Ati/N, independent of the order k of represen- 
tation. We observe that v s is double of that appearing in the homogeneous model (k = 1) 
H| 14] . This fact can be easily interpreted, by noticing that the elementary translation 



"cell" of the alternating mixed models is double of that of homogeneous system. Indeed, 
considering this discussion and the previous results of ref. [13|, [14|] (assuming independence 



of k) will lead us to guess that v s = An/N, with no need of an explicit computation 0. 

In the next section we are going to compute the conformal anomaly which defines the 
class of universality of these conformal mixed GL(N) models. 

3 The thermodynamics of the mixed GL(N) model 

In order to discuss the thermodynamic properties we adopt the thermodynamic Bethe 
ansatz approach originally proposed by Yang and Yang [15|| . This method is based on the 
minimization of the free energy and takes advantage of the integrability through the Bethe 
ansatz equations. The first step is to notice that the Bethe ansatz equations ([7]) admit 



the same string hypothesis used previously by Takahashi |16| in the isotropic Heisenberg 
model. This observation allows us to conclude that , in the L — > oo limit, the parameters 
X r j organized in strings of the type described in Eq([L6|). Substituting Eq.(|l^) in Eq.(|7|) 
and taking the thermodynamic limit, we are able to obtain the following infinity set of 
coupled integral equations for the densities cx^(A) (<7^(A))of particles (holes) 



8 1 r N-l oo 

*»< A ) = j£ M^) + <W(A)] - E E(4* * B r y * a] )(A) (17) 

r'=li =1 



2 In fact, for an alternating model with periodicity / ( made by a collection of "spins" operators at 
different order of representation), we should have v s = 2lir/N . 
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where n indicates the length of the n-string, and (/ * g)(x) denotes the convolution 
2^ J^oo f( x ~ y)9(y)dy- The functions A nj (A), B rr >(\) and </> n j(A) are easily represented 
in terms of their Fourier transforms. Defining the Fourier component of a given function 



f(x) by f(u>) = 5- dxe lulx f(x), we have the following expressions for A n j(u>), B r 
and cj) n j(uj) 



A nJ (uj) = coth(|u;|/2) _ e -(r i+J -)|^|/ 2 

B r y(u) = 6 r y - p(u)l r y 

4>n,j/2{u) = Anj(u)p(u)) 

and I > is the incident matrix of the A N _i Lie algebra. 



A") 

(18) 
(19) 
(20) 



where p{uS) 



2cosh(u/2) 

The second step is to encode the temperature T via minimization of the free en- 



ergy F sym = E sym — TS. By using a standard procedure |u| [16], [u], [18] the energy 
E sym and the entropy S can be written in terms of the densities of particles (<r£(A)) and 
holes (cr^(A)). After the minimization, SF sym = 0, we get the following thermodynamic 
Bethe ansatz (TBA) equations 



N-l 



<(A) 



K r (X)(5 n>1 + 6 Ujk ) <P r ,r' * 



ln(l + e e " /T ) 



N-l 



T <Prr' * M 1 + e £ "+ l/T ) + ln(l + e<-i /T ) 



(A) 
(A) 



(21) 



where <7^(A)/er^(A) = e e ™^' T , 6q(A) = 0, and the Fourier component of the functions 
K r {uj), (p rr >(uj) and rr >(\) are given by 



K r {uj) 

<PrS M 



pHb-1(u) 

5„' — B , , 



p{u)B-\ 

r.r 



(22) 
(23) 



Finally, the equilibruim free energy F sym is given by 



^sym 



T roo N-l 

/L = e sym - — d\Y] KJX) ln(l + e^' T ) + ln(l + eW T ) 

4:71 J-OO ^ L 



(24) 
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One possible way to calculate the central charge of a conformally invariant system is 
by analysing the low temperature behaviour of the respective free energy. The universal 
behaviour of the free ene rgy is given by 

F/L = eoo - ^ (25) 

It turns out that Eqs(pT|, |24] ) allow us to make an exact calculation of such low temper- 
ature behaviour. We first define the shift A — > A — ^- ln(^j-), taking the derivative in A of 
Eq. (|2T|) and after some few standard manipulations |17j the T — > limit can be expressed 
in terms of the Dilogarithm functions L(x) by 



F sym /L = e 



sym 



NT 2 
^2A 



N-l fc-2 



(tf-i)*-E E% 



sin[(A; — m — 1)9] sm[m9] 



r=1 m=1 l sin[(m + r)9] sin[(r + k — m — 1)9] ' 
where 9 = ^ and L{x) = _^ jj ^[MJ + 



(26) 



Using some identities for the sum of the Dilogarithm function proved in ref. [^T], we 
finally have 

F sym /L _ eSym T 2 (N-l)((N + 2)k-2) 

1 600 24 N + k-1 [n 

Comparing Eqs. ([25], |27| ), we find that the central charge is c = (iV — l)((iV + 2)fc — 
2)/(N+k—l). Remarkably enough, this conformal anomaly can be decomposed in terms of 
the central charges of two SU(N) WZWN models with topological charge k — 1 (c = N—l) 
and k — k — 1 (c = (iV 2 — l)(fc — 1)/(N + k — 1)). This result generalizes similar 
decomposition mentioned by the authors || for the SU(2) mixed Heisenberg model. In 
order to give an extra support for this value of the central charge, in appendix B we present 
some numerical results for the finite size effects of the ground state energy. 



4 Discussions on possible generalizations 

It is almost evident that all of our discussion of section 2 can be generalized to an arbitrary 
representation of order k in the auxiliary space of states. The main technical difficulty 
is the explicit construction of the non-isomorphic vertex operator •(//). The solution 
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of this problem has already been considered in ref. || for arbitrary finite representations 
of GL(2) Lie algebra. The vertex .(//) is expressed as a linear combination of certain 
projectors defined on the subspaces of the Klebsch-Gordon decomposition of GL{2)k ® 
GL(2) k >. The transfer matrix T kk >(fi) is then defined by 

= Tr v{k) [Rl^)R k k , ^{v) ■ ■■R k k>2 ( f i)R k klil ( f i)} (28) 

The eigenvectors and the eigenvalues of the associated one dimensional Hamiltonian 
H k k > can be determined by using the following strategy. We first define an auxiliary 
transfer matrix T^f(/j,) which commutes with T kk '(fi) as 

T fc 7(//) = Tr vW [RlMR\, ^) ■ ..Rl^Rl, ^)} (29) 

This definition has the advantage of reducing the auxiliary space to its fundamental 
representation, hence the standard QIS formalism can be applied. On the other hand, 
the eigenvalues of T k k >(n) can be related to those of T""f(/i) by certain commutators of 
T fcfc /(/i) and the usual quantum inverse scattering B(n) operator |T7| . The eigenvalues of 
the H k k < are parametrized by the Bethe ansatz equation 

( \-ik/2 \ L/2 ( X s -ik'/2 \ L/2 = _*L Xj-^-i 

\\j + ik/2 ) \Xj + ik'/2) fJi Xj - Xi + i 1 1 

where /ij = iXj — 1/2 and the eigenvalues E k fc / are determined by 

M k 

^ = ~§(Ai) 2 + (fc/2) 3 (31) 

Analogously to what we have discussed in section 2, similar approach works for T k > fe (/i) 
and the conformally invariant quantum Hamiltonian can be defined through the product 
^WCaO^V *;(/■*)■ Considering the results of section 2 and those of ref. [f| on the Yang- 
Baxter solutions for the GL(N) group, it seems plausible that similar conclusions reached 
for N = 2 can be extended for an arbitrary value of N. Comparing the left hand side of 
Eqs(^, |30|) in the case of iV = 2, we observe that a certain factor 1/2 has been replaced by 
k/2. This leads us to conjecture that similar mechanism should work for a general mixed 
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GL(N) system. Taking into account this fact and writing Eq.(^) in a more convenient 
way, we conjecture that the form of a GL(N)k ® GL(N) k > Bethe ansatz equation is 

f \)- i 5 7 , l k/2 \ L ' 2 f y. - iS^k' /2 \ L/2 = _ N n 1A ^ \j -iC r y/2 

\X5 + iSr,ik/2) {\r + i5 rjl k' /2 ) -X\ +iC r y/2 [ ' 

where C rr > is the An-i Cartan matrix and the eigenvalues E sym of the conformally invari- 
ant Hamiltonian are given by 

F sym \^ * \ ^ Z (<¥>\ 

^(Aj) 2 +(w h ( A )) 2 + (*v 2 ) a 

It is not a surprise that the structure of the Bethe ansatz equation is closely related to 
the An-i Lie algebra. In the case of homogeneous vertex models (k = k), the 

authors of ref. |S| have conjectured that the same structure will remain for all semi- 
simple A,D,E Lie algebras^. In fact in ref. |23] this conjecture has been verified by an 
explicit computation in the case of D n Lie algebra. Based on these observations, let us 
assume that the same conjecture can be extented to the case of non- homogeneous (k®k') 
models. It is not difficult to verify that the associated TBA equations are similar to 
the system of Eqs. (|2ll) . We just have to replace l rr > and iV in Eqs. fl2lD by the incident 
matrix and the rank of the corresponding A,D,E Lie algebra. Interesting enough, these 
equations can be cast in a rather useful form which will be helpful in the analysis of the 
low temperature . Defining the function Y£(\) = e _e ™^' T , making few manipulations 
in the Fourier transform of Eq.([2lD and Fourier transforming back we find the following 
expression 

-i 



r eG 



Y;(\+i/2)Y:(x-i/2) n i+^ r (A) rr n V + ^/y;w 



ln,m _ 27r<5(A)<5 r ,i(5„, fc +5 .,)/T 



(34) 

where r is an index characterizing the nodes of the Dynkin diagram of the G = A, D, E 
Lie algebra and j is a similar (unrestricted) index for the A^ Lie algebra. 



3 This fact is related to the idea that the classification of the solutions of the Yang-Baxter equations 
is somewhat connected to the classification of the Lie algebras and their automorphisms. 
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Eq. (|34D defines a set of functional hierarchy relations for the functions Y*(X). The 
possibility of constructing such functional relations from the TBA equations has been 
first noted by Al. Zamolodchikov |?5] in the case of the diagonal system of scattering 
S'-matrices. It also appears that certain functional hierarchies play the keystone in the 
computation of critical dimensions in the integrable lattice models [26fl. In our case they 
encode all necessary information in order to obtain the low temperature behaviour of the 
free energy. Proceeding as in sect. 3 , we can show that the T — > limit of the free energy 
assumes the following form 



F sym /L = e 



sym 
oo 



T 2 h G 
24tt 



k'ra-Y, E £(1/1 + #(*))-£ E 

reGjeA k r£Gj£A k ,_ k 

(35) 



where r G and h G are the rank and the dual Coxeter number of the Lie algebra G, and 



B -i 

e sym = _ joo^ du 4cos ^l /2) [<fa |fc/2 (u) +</> k ' >k '/ 2 M +2^, fc/2 H] ■ The constants j/J(m) satistfy 



oo 



the equation 



n 

r'eG 



1 + j/J (m) 



n [l+iMw 



where Z rr ' (i^-) is the incident matrix of the Lie algebra G (A m ). 

Remarkably enough, the sum of Dilogarithm function |2T]. [22j appearing in Eq. 
been conjectured to have the expression 

roTn(m — 1) 



E E L(l/l + fJm)) 



ha + m 



(36) 



has 



(37) 



reGjeA m ,L G 

Although the proof of last identity was essentially given for the A and D Lie algebra, it 

is quite direct to verify it by numerically solving Eq. ([36D for several small values of m and 

ho |p2[ |28|| . Using this Dilogarithm sum and taking into account that now v s = An /he, 

we finally obtain the following central charge 

r G k(h G + l) + r G (k' -k)(h G + l) 
h G + k h G + k' — k 

where we have already decomposed the result in terms of the central charge of two G 
invariant WZWN models with topological charges k = k and k = k — k, respectively. 
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In particular for k = k we recover the known conjecture that rational isomorphic vertex 
models are in the class of universality of WZWN field theories [|, [2Cj, |28| . 



5 Conclusion 

In this paper we have discussed the critical behaviour of conformally invariant mixed 
GL(N) spin chains. The central charge of an alternating model with "spins" operators 
at order k (k') of representation acting on even (odd) sites can be decomposed in terms 
of WZWN field theories with topological charge k and k' — k (k' > k). We have also 
considered the generalization of this result to other symmetric representation of the A, D, E 
Lie algebras. 

Another possible extension of the results of this paper is to consider a collection of ver- 
tex operators at different representation distributed on a line of size L and with periodicity 
I. The associated rotational symmetric transfer matrix can be defined as 

-'Zl'jJin = I[Tr^\R^M--.Rti L _U^ (39) 
t=i 

Following our considerations of sections 3 and 4, the basic change in the TBA equations 
(34) is that the right hand side of Eqs. ([54]) is replaced by 2Tc5(X)S rj i Ya=i S n ,kjT. For 
instance, taking the following ordering k\ < k 2 < ■ ■ ■ < ki, the central charge of the one 
dimensional Hamiltonian associated to the system (^) will be 

i-i 

c = -A*); k = (40) 

i=0 

where c(k) is the central charge of the A, D, E WZWN model with topological charge k. 

Finally, it should be interesting to study the full operator content of these models 
and in particular to understand the decomposition of Eq.(pE0"D in terms of the bosonic and 
parafermionic fields of the WZWN theories. We notice that for the sequence k i+ i — ki — 1 
and k\ = 1, the central charge is / time the rank of G and hopefully the operator content 
will be determined by Ire coupled bosonic fields. 
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Appendix A 



Following the basic steps of the quantum inverse scattering method we propose a set 
of eigenstates \ip > defined by [@, §, |ll] 



|V ^04, ■ ■ A*i) = F jv .. jMl BK(ti) ■ • -5^(^)10 > (A.l) 

where |0 > is the pseudo vacuum (see Eq.(5)). The next step, motived by the properties 
of Ti ik (/i)\0 >, is to decompose the monodromy matrix as 

n,M=( AM B ' M ) (A.2) 

^ C>(jm) D»(ji) ) 

where is a row(column) index, i, j — 1, - • • , N — 1. 

It follows from the identity Ti^/i) = Tr^i) [^^(/i)] that 

T^IV >= [ACaO + X^O*)]^ >= A(M/4})I^ > (A.3) 

i 

In order to find the eigenstates \ip > and the eigenvalues A(/i, {//[}) we also need the 
commutation relation between the B l (fi) and D n (fi) operators. These relations 

follows from the relation 

R\(fx - n)r 1:k (n') ® r Xifc (/i) = r ljfc (^) <g> Ti ifc (//).&[(// - /i) (A.4) 



where = VR\(fi). Using Eqs.(pj[A2|, gg) we have 



1 1 

A(ai')S*(aO = ± ^(/i)A(/i') + F^'jA^) 

^tX ^(X ^(X ^Li 
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where [R\(fi — are ^ e GL(N — 1) matrix elements of the matrix defined in Eq.(||). 

From Eqs.( [01 , |X^ , |70| ) follows that 



Ml u 1 - a + 1 
A{^> = II i o(^)lV»>+U.T. 

f=i f4 - t* 

M 1 i . . 

1 , 



^iV* > = (II Tl)^H--C Fll ''' tMlBl1 ^ ■ ■ ■ B ' Ml (^mO^)I^ > +U.T. 



8=1 



// // ' ' M/ 



(A.6) 



where a(/i) = (/i + l) L / 2 (/i + |(Jfe + 1) L / 2 , d(/x) = /i L / 2 (/i - \{k - l)) L / 2 and t%?£[ are 
matrix elements of the following operator 

t (2) (/i; {//)}) = tr[R{ jMl (n - ^[n-i] ■ • ■ #i,iO" " /4)[iv-i]] (A.7) 

and F h "" 1 m^ are eigenvector's component of t^(fi) with eigenvalues A( 2 )(/i, {/i*}). 

The U.T. symbol stands for "unwanted terms" which appears due to the interchange 
of the arguments /x and \i in the relation ( |A.5|) . When these terms are null, \ip > becomes 
an eigenstate of T 1)fc (/i) and as a consequence we obtain a restriction to the rapidities 
H (the Bethe ansatz equation). Finally, Eq.( A.6 ) is solved by introducing in each step 
i = 2, • • • , N — 1 a new matrix t^(/i) acting on sites, analogously to that of Eq. flA.7p . 



The final result for the eigenvalues A(p) of Ti t k(p) and A^/i, 1 }) of t^(fx) are 

A/ 1 1 _ , I Af 1 i 

A(ai) = a(/i)II - / + <*(aOII rA^^^f}) (A.8) 

AT"- 1 AT r _ , i Af- 1 A/ r i 

A (r V; i^- 1 }) = n - Mr 1 + 1) n ,/, + n ^ - ^r 1 ) n -4-7 a^v; k» 

i=l j=l r i=l i=l Z 1 Z 1 * 

Eq.(7) is then obtained by imposing the zero residue condition in Eq.(A.8). 



Appendix B 

The critical behaviour of a conformally invariant theory can be determined by studying 



the consequences of the finite size L effects for the eigenspectrum [27]. For example, the 
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central charge is related to the ground state energy E sym (L) by JT9|, ^0 



7T1) C 

E sym (L)/L = eZ m - -j f- (B.l) 

The central charge c can be numerically calculated by extrapolating the sequence 

c(L) = —[E sym (L)/ L - eZ" 1 }^ (B.2) 

In table l(a,b) we present our estimatives for the sequence (B.2) in the case of iV = 2, 3 
( The N = 2 data has been already presented by us in ref.||.) and k = 3. We notice that 
these numerical results are in concordance with the TBA analysis of section 3. (Eq.(27)). 
In our numerical analysis we have also observed that the case k = 2 is rather special. The 
string hypothes is (X r j = £J ± i/2) is almost exact for large enough L, presenting a very 



unusual small correction |j, In this case we can use the analytical method of ref. [29| 
and conclude that the central charge is c = 2(A^ — 1) (in agreement with Eq.(27)). 



Table l(a,b) The estimatives of the central charge of Eq.(B.2) for k = 3, (a) N = 2 
and (b) N = 3. 



(a) 



(b) 



L 


N = 2,k = 3 


8 


2.839 364 


16 


2.602 543 


24 


2.556 301 


32 


2.538 530 


40 


2.529 503 


48 


2.524 142 


Extrapolated 


2.500(6) 



L 


N = 3,k = 3 


12 


5.614 431 


24 


5.336 131 


36 


5.277 662 


48 


5.254 475 


60 


5.242 437 


72 


5.235 167 


Extrapolated 


5.206(1) 
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